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We consider a self-consistent system of Bianchi type-I (BI) gravitational field and a binary 
mixture of perfect fluid and dark energy. The perfect fluid is taken to be the one obeying the 
lO ! usual equation of state, i.e., p = l^e, with G [0, 1] whereas, the dark energy is considered 

to be obeying a quintessence-like equation of state. Exact solutions to the coiTcsponding 
^ ■ Einstein equations are obtained. The model in consideration gives rise to a Universe which 

■ is spatially finite. Depending on the choice of problem parameters the Universe is either 

close with a space-time singularity, or an open one which is oscillatory, regular and infinite 
in time. 
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in ; I. INTRODUCTION 



The discovery that the expansion of the Universe is accelerating IjJ has promoted the search 
for new types of matter that can behave like a cosmological constant 0111] by combining positive 



O 

$—1 . 

bJQ. energy density and negative pressure. This type of matter is often called quintessence. Zlatev et 
at. 1 4] showed that "tracker field", a form of qiuntessence, may explain the coincidence, adding 
■ new motivation for the quintessence scenario. 
5_j . An alternative model for the dark energy density was used by Kamenshchik et at. where the 

authors suggested the use of some perfect fluid but obeying "exotic" equation of state. This type 
of matter is known as Chaplygin gas. In doing so the authors considered mainly a spatially flat, 
homogeneous and isotropic Universe described by a Friedmann-Robertson-Walker (FRW) metric. 

The theoretical arguments and recent experimental data, which support the existence of an 
anisotropic phase that approaches an isotropic one, lead to consider the models of Universe with 
anisotropic back-ground. Since the modem-day Universe is almost isotropic at large, its simplicity 
and evolution into a FRW Universe makes the BI Universe a prime candidate for studying the 
possible effects of an anisotropy in the early Universe on present-day observations. In a number 
of papers, e.g., |0,|7|l, we have studied the role of a nonlinear spinor and/or a scalar fields in the 
formation of an anisotropic Universe free from initial singularity. It was shown that for a suitable 
choice of nonlinearity and the sign of A term the model in question allows regular solutions and the 
Universe becomes isotropic in the process of evolution. Recently Khalatnikov et al. 1 8] studied the 
Einstein equations for a BI Universe in the presence of dust, stiff matter and cosmological constant. 
In a recent paper |@] the author studied a self-consistent system of Bianchi type-I (BI) gravitational 
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field and a binary mixture of perfect fluid and dark energy given by a cosmological constant. The 
perfect fluid in that paper was chosen to be the one obeying either the usual equation of state, i.e., 
p = C,e, with C, E [0, 1] or a van der Waals equation of state. That paper was followed by another 
where we studied the evolution of an initially anisotropic Universe given by a BI spacetime and a 
bimnary mixture of a perfect fluid obeying the equation of state p = L,£ and a dark energy given 
by either a quintessence or a Chaplygin gas HlOl] . It should be mentioned that the inclusion of 
dark energy does not eliminate initial singularity of the model and the space-time in those cases 
is ever-expanding. In order to obtain a singularity-free Universe in the present paper we introduce 
a modified version of quintessence-like dark energy, which at the same time is able to explain the 
accelerated expansion. 

II. BASIC EQUATIONS 

The gravitational field in our case is given by a Bianchi type I (BI) metric in the form 

ds^ = dt^ - a^dx^ - b^dy^ - c^dz^, (2.1) 

with the metric functions a, b, c being the functions of time t only. 

The Einstein field equations for the BI space-time we write in the form 
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Here k is the Einstein gravitational constant and over-dot means differentiation with respect to t. 
The energy-momentum tensor of the source is given by 

T; = {e + p)u^u''-p5;, (2.3) 

where is the flow vector satisfying 

g^yu^u'' = \. (2.4) 

Here e is the total energy density of a perfect fluid and/or dark energy density, while p is the 
corresponding pressure, p and £ are related by an equation of state which will be studied below in 
detail. In a co-moving system of coordinates from ( 12.31 ) one finds 

= £, r/ = Ti = = -p. (2.5) 

In view of (12.51) from (12.21) one immediately obtains |0] 

/I. 
(2.6a) 

/I. 
■^], (2.6b) 

c(t) = D3Ti/3exp[X3 J (2.6c) 
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Here D,- and X,- are some arbitrary constants obeying 

DiD2Z)3 = l, Xi+X2+X3=0, 

and T is a function of t defined to be 

T = abc. (2.7) 
Summation of (I2.2al) . (I2.2bl) . (I2.2cl) and 3 times (I2.2dl) gives the equation for T : 

I-Yi^-P)^ (2-8) 
whereas, from the conservation law for the energy-momentum tensor for e we find 

e = ~{e + p). (2.9) 



After a little manipulations from (12.81) and ( 12.91) we find 

t = ±^/Cl+3K£T^, (2.10) 
with Ci being an integration constant. Eq. ( 12.91) can be rewritten in the form 



T 



(e + p) x' 



(2.11) 



Taking into account that the pressure and the energy density obey a equation of state of type 
p = /(e), we conclude that e and p, hence the right hand side of the Eq. (12.81) is a function of T 
only, i.e., 

z = ^{£-p)t = J^{t). (2.12) 

From the mechanical point of view Eq. (12.121) can be interpreted as an equation of motion of a 
single particle with unit mass under the force ^{t). Then the following first integral exists llllh : 



t= ^2[^-^(t)]. (2.13) 

Here S' can be viewed as energy and ^(t) is the potential of the force Comparing the Eqs. 
(irm and dTH) one finds = Ci /2 and 

^{t)^~k£t^. (2.14) 

Finally, rearranging (12.101) . we write the solution to the Eq. (12.81) in quadrature: 

dz 

, =t + tQ, (2.15) 

VCi+3K:eT2 

where the integration constant Iq can be taken to be zero, since it only gives a shift in time. 

In what follows we study the Eqs. (12.81) and ( 12.91 ) for perfect fluid and/or dark energy for 
different equations of state obeyed by the source fields. 
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III. UNIVERSE AS A BINARY MIXTURE OF PERFECT FLUID AND DARK ENERGY 

In this section we thoroughly study the evolution of the BI Universe filled with perfect fluid 
and dark energy in details. Taking into account that the energy density (£) and pressure (p) in this 
case comprise those of perfect fluid and dark energy, i.e., 

e = £pf + £DE, P = Pp{ + PDE 

the energy momentum tensor can be decomposed as 

= {£DE + £pf + PDE + Ppf)u^U^ - (pDE+Ppf)S^- (3.1) 

In the above equation Cde is the dark energy density, />de its pressure. We also use the notations 
Cpf and ppf to denote the energy density and the pressure of the perfect fluid, respectively. 

In a comoving frame the conservation law of the energy momentum tensor leads to the balance 
equation for the energy density 

£de + £pf= — [£de + £pf + Pde + Ppf) ■ (3.2) 

The dark energy is supposed to interact with itself only and it is minimally coupled to the gravi- 
tational field. As a result the evolution equation for the energy density decouples from that of the 
perfect fluid, and from Eq. (13.21) we obtain two balance equations 

t / ^ 

Cde + - (cde + Pde) = 0, (3.3a) 

T 

epf+-(epf + ;?pf) = 0. (3.3b) 

'V 

A. Equations of state 

In order to complete the system of equations we need to specify two equations of state for ppf 
and Pde- 



1. perfect fluid 

There are some equation of states that are commonly used that, although not widely applicable, 
are obtained as a result of approximate estimates for particular fluid. The barotropic equation of 
state 

Ppf =C Cpf, (3-4) 
is often assumed. Here ^ is a constant and lies in the interval ^ G [0, 1]. Note that < ^ < 1 is 
necessary for the existence of local mechanical stability and for the speed of sound in the fluid to 
be no greater than the speed of light. Depending on its numerical value, ^ describes the following 
types of Universes lfl2ll 

^ = 0, (dust Universe), (3.5a) 

^ = 1/3, (radiation Universe), (3.5b) 

C G (1/3,1), (hard Universes), (3.5c) 

^ = 1, (Zel'dovich Universe or stiff matter) . (3.5d) 
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In view of the Eq. (13.41) from (I3.3bl) one easily finds laws of change of energy density and 
pressure of a perfect fluid with the expansion of the Universe : 

£,f = eo/T^'+^\ Ppf = eoC/r^'+^\ (3.6) 
where Cq is the integration constants. Li absence of the dark energy one immediately finds 

T = a2/(i+o, (3.7) 

with C being some integration constant. As one sees from (12.61) . in absence of a dark energy, for 
^ < 1 the initially anisotropic Universe eventually evolves into an isotropic FRW one, whereas, 
for C = 1' i-e., in case of stiff matter the isotropization does not take place. 



2. Dark energy 

It was mentioned earlier that the dark energy can be given by a A term, a quintessence or a 
Chaplygin gas. It was also mentioned that the quintessence was constructed by combining positive 
energy density and negative pressure and obeys the equation of state 

Pq = WqCq, (3.8) 

where the constant Wq varies between —1 and zero, i.e., Wq G [— 1, 0]. 

Here we introduce a modified model of quintessence when the dark energy and the correspond- 
ing pressure obeys the following equation of state: 

PDE = -w(eDE - Ccr) , (3.9) 

where the constant w G [0, 1). Here Cd- some critical energy density. Setting £cr = one obtains 
ordinary quintessence. It is well known that as the Universe expands the (dark) energy density 
decreases. As a result, being a linear negative function of energy density, the corresponding pres- 
sure begins to increase. In case of an ordinary quintessence the pressure is always negative, but 
for a modified quintessence as soon as Cq becomes less than the critical one, the pressure becomes 
positive. In Fig. [l]we illustrate the evolution of pressure corresponding to ordinary and modified 
quintessence with the expansion of the Universe. 



(3.10) 



In account of (13.91 ) from (I3.3al) one finds the following relation between £de and T 

1 



Cde = 

I — w 

with El being some integration constant. 



t1 



B. Exact and numerical solutions 

As soon as the right hand side of the Eq. (12.121) is defined, we can study this equation in details. 
First we write the solution to the equation in question in quadrature which will be followed by 
numerical results. 

Inserting Cpf and Eq into (12.151) one now finds 

/ , =t + to. (3.11) 

^ /Ci +3K(eoT(i-^) + [ei/(l -w)]t(^+^) - [we^Jil -w)]t^) 
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FIG. 1: Evolution of the pressure with the expansion of the BI Universe when it is filled with perfect fluid 
and dark energy obeying ordinary and modified quintessence given by ( I3.8t and ( I3.9I ). respectively. Here 
the letters "o" and "m" stand for ordinary and modified, respectively. 



Here to is a constant of integration that can be taken to be trivial. As one sees, the positivity 
of the radical imposes some restriction on the maximum value of T, i.e., in this case the model 
allows oscillatory mode of expansion. It means the dark energy initially acts as a repulsive force 
resulting in accelerated expansion of the Universe. But as soon as Cde becomes less than the 
corresponding pressure changes its direction and as a result the Universe begins to contract. In 
order to give a complete picture, beside the system with a modified quintessence we study the 
system with an ordinary one as well. It should be noted that the critical energy density Ede should 
be very small. Here for simplicity we set = 0.01. 




FIG. 2: View of the potential when the BI Uni- 
verse is filled with a binary mixture of perfect 
fluid and a quintessence given by (13.81) . 



FIG. 3: View of the potential when the BI Uni- 
verse is filled with a binary mixture of perfect 
fluid and a modified quintessence-like dark 
energy given by (13.91) . 



In Figs. El and |3l potentials corresponding to an ordinary and modified quintessence are given. 
As one sees, the usual quintessence does not allow oscillatory mode of expansion. Universe in this 
case expands endlessly. A corresponding expansion of the Universe is given in Fig. |4| It should be 
emphasized that the initially anisotropic Universe in this case evolves into an isotropic FRW one. 



Anisotropic cosmological models with perfect fluid and dark energy revisited [February 7, 2008J7 



As oppose to the ordinary quintessence a modified quintessence imposes some restriction on the 
maximum value of T. As a result the Universe initially expands, but after reaching some maximum 
it begins to contract. Depending on the choice of the constant, which can be viewed as an energy 
level, it either shrinks into a point [cf. Fig. 13 thus giving rise to space-time singularity, or begins 
to expand again after reaching some non-zero minimum, i.e., the Universe in this case experiences 
the oscillatory mode of expansion [cf. Fig. Note that to each S' corresponds a particular pair 
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FIG. 4: Evolution of the BI Universe is filled 
with a binary mixture of perfect fluid and a 
quintessence. For simplicity, as a perfect fluid 
we consider only radiation. 



FIG. 5: Evolution of the BI Universe with a 
modified quintessence. Choosing (f > (here 
we set = 1 one obtains nonperiodic picture 
of evolution. As one sees, different choice 
of gives rise to different amplitude, but the 
overall character of the solution remains unal- 
tered. 
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FIG. 6: For < BI Universe filled with a bi- 
nary mixture of perfect fluid and quintessence- 
like dark energy admits oscillatory mode of 
expansion. Here we set = —1000 with 
To = 50. 



FIG. 7: View of energy density and pressure 
when BI Universe experiences oscillation as 
in Fig. |6l 
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In Fig. |7]we illustrate the evolution of energy density and pressure when the Universe is filled 
with a binary mixture of perfect fluid given by a radiation and a modified quintessence. As one 
sees, the initial density is large enough and the pressure is initially negative, with the energy 
density being less than £cr, pressure becomes positive and the Universe begins to contract. As a 
result the energy density begins to increase and pressure again becomes negative. This results in 
the oscillatory mode of expansion. 



70000 




t 



FIG. 8: View of the acceleration for different source fields. Here "rad", "quint" and "quint-m" stand for 
radiation, a mixture of radiation and an ordinary quintessence and a mixture of radiation and modified 
quintessence, respectively. 

Finally in Fig. [5] we plot the graphic of acceleration verses time. As one sees, if the Universe 
is filled with perfect fluid only, the process of evolution decelerates with time, while the inclusion 
of an ordinary quintessence results in an accelerated evolution. If the dark energy is given by a 
modified quintessence the accelerated expansion is followed by a decelerated one. 

IV. CONCLUSION 

A self-consistent system of BI gravitational field filled with a perfect fluid and a dark energy 
given by a modified version of quintessence has been considered. The exact solutions to the 
corresponding field equations are obtained. The inclusion of the dark energy into the system gives 
rise to an accelerated expansion of the model. As a result the initial anisotropy of the model 
quickly dies away. The modification of the quintessence results in appearing some upper bound of 
the volume scale, i.e., in this case the Universe is spatially finite. Depending on the choice of the 
integration constant S', which can be viewed as an energy level, the Universe is either close with a 
space-time singularity, or an open one which is oscillatory, regular and infinite in time. 
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